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ABSTRACT 



Asteroseismology of solar-type stars has entered a new era of large surveys 
with the success of the NASA Kepler mission, which is providing exquisite data 
on oscillations of stars across the Hertzprung-Russell (HR) diagram. From the 
time-series photometry, the two seismic parameters that can be most readily 
extracted are the large frequency separation (Az/) and the frequency of maximum 
oscillation power (z^max)- After the survey phase, these quantities are available 
for hundreds of solar-type stars. By scaling from solar values, we use these two 
asteroseismic observables to identify for the first time an evolutionary sequence 
of I-M0 field stars, without the need for further information from stellar models. 
Comparison of our determinations with the few available spectroscopic results 
shows an excellent level of agreement. We discuss the potential of the method for 
differential analysis throughout the main-sequence evolution, and the possibility 
of detecting twins of very well-known stars. 

Subject headings: asteroseismology — stars: evolution — stars: oscillations 



Introduction 



The high-precis ion photometric observations obtained by the NASA Kepler Mission 
( Borucki et al.ll2009l ) have led to a dram atic increase in the number of main-sequen ce and 
subgiant stars with detected oscillations (IChaplin et al.ll2010l ICilliland et al.ll2010al ). Data 
on hundreds of these stars are now available to test our knowledge of solar-like oscillations 
and their dependence on stellar parameters. Moreover, extremely precise information can be 
extracted on fu ndamental stellar prop erties and compared to population synthesis models of 
our galaxy (e.g. IChaplin et al.ll2011al ). 



The large ensemble of stars with detected oscillations allows us to select cohorts of 
targets sharing one or more common properties. By performing comparative studies on 
these stars, it should be possible to suppress the dependence of the modeling results on the 
shared characteristic, opening an exciting possibility for further constraining the internal 
physical processes. This is the basis for differential asteroseismology: a detailed seismic 
comparison of the inner structures of stars showing some similar property. 

To exploit this technique, one must be able to identify stars with common characteristics. 
Comparing the quantities extracted from conventional observations such as spectroscopy is 
not sufficient to fully constrain stellar properties, and evolutionary models must be employed 
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to determine e.g. masses and radii. Asteroseismology allows us to go one step further by 
providing a tool to find pairs or groups of stars with similar characteristics. 

In this Letter, we use global seismic properties of solar-type targets observed by Kepler 
to construct an observational sequence of field stars with masses similar to that of the Sun, 
in a model-independent way. 



2. Extracting mass and radius from asteroseismic data 



Two asteroseismic parameters can be readily extracted from the p-mode oscillation spec- 
trum without the need for individual frequency determinations. One is the large frequency 
separation, Au, which denotes the frequency difference between modes of the same degree 
and consecutive radial order. The other is the frequency of maximum oscillation power, 
^'max- It has been shown that, to ver y good approximation, Au scales as the square root 
of the mean density (e.g. lUlrichI 119861 ). while z/^ax is related to the dyn amical timescale of 



the atmosphere, as repres ented by the acoustic cut-off frequency (e.g. iBrown et al.l Il991 
Kjeldsen fc Bedding! Il995l ). These two quantities are tightly correlated over a wide range 
of values, arid follow scaling relatioiis from the accuratel y known solar parameters (e.g. 
Hekker et al.ll2009l . IStello et al.ll2009al . iHekker et al.ll2011al ). These scaling relations can be 
written as 

Au 
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where T^s is the effective temperature of the star, and A vj^ = 135.1 /xHz, t^max,© = 3150 /iHz, 
^eff,0 = 5777 K are the observed values in the Sun (e.g. IChaplin et al.ll2011b[ ). From these 
two equations it is evident that, provided we have a measurement of Tes, the global seismic 
observables give a determination of stellar mass and radius for a given star that is independent 
of evolutionary models. This is the so-called direct method, and it has been used to constrain 
the fundamental properties of red giant stars and distinguish different populations of stars 



Miglio et al. 



20091 . iKallineer et al.ll2010l . iHuber et al.ll2010l . iHekker et al.ll2011b[ ). It can 



(e.g. 

now be applied to the hundreds of solar-type stars for which Kepler has detected oscillations 
Jchaplin et allboilal ). 
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Fig. 1. — T'eff-t'max diagram for the complete sample of targets with detected oscillations. 
Evolutionary tracks are also plotted: at solar ([Fe/H]=0.0) metallicity for 0.85 Mq (dash- 
dotted line), 1.0 Mq (solid line) and 1.15 M0 (dashed line), and at sub-solar metallicity 
([Fe/H]=-0.5) for 1.0 Mq (dotted line). The Sun is marked with its usual symbol. 
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3. Data analysis and stellar properties determination 

We examined asteroseismic results for 2000 solar-type stars observed by Kepler during 



Verner et al. 


2011 


). Each s' 


.ar was 


sampling; see 
described by 


Gilliland et al. 


2010b) 


Garcia et al. 


( 


2011 





2010b[) . Time series were prepared 



on the raw light curves. The prepared light curves were analyzed by different teams, who 
detected and extracted the basic properties of the solar-li ke oscillations. The methods of ex- 
traction and descriptions of each pipel i ne may be found in IChapl in et al.l (120101). iHuber et al.l 



feoioh 



and 



Karoff et al.l (1201 ll ). 



(jgOOgh . lMosser fc AppourchauxlJ2009h . lHekker eTaP J201oh . lMathur et al.lfcoioh . lCampante et al 



Effect ive temperatures have b een derived via the InfraRed Flux Method (IRFM) de- 
scribed by ICasagrande et al.l (l2010l ). using the available Sloan and 2MASS photometry com- 
piled for the targets in the Kepler Input Catalog (KIC; iBrown et al.ll201l[ ). 



Figure H] shows a Teff-i^max diagram constructed from the consolidated results of different 
pip elines for all the ta rgets with oscillations measured by at least two of them (as described 
by IVerner et al.ll201ll ). Stars with z/max below 350 /iHz are beyond the subgiant phase of 
evolution and have been omitted from our sample. We have overplotted evolutionary tracks 
for different masses and metallicitie s calculated with the GARching STe llar Evolution Code 
fGARSTEC. IWeiss fc Schkttllboosh . scaled from the lCrevesse fc Sauvallf 
considering a helium enrichment ratio of AY/AZ = 1.8 ( ICasagrande et al. 



1998) solar mixture 



20071 ) and not 



including diffusion of helium and heavy elements (which is why the solar metallicity 1-Mq 
track does not overlap the position of the Sun). 

It is already clear in Fig. [1] that a number of oscillating stars are located on or very close 
to the evolutionary tracks and thus are potentially 1-Mq stars. Interestingly, the stars are 
concentrated in two main regions, with an apparent gap between them (around T^s ~5300 
K and I'raa.x ~600 /iHz). This feature suggests a rapid evolutionary stage where few stars are 
found, which could be identified as the transition between the main-sequence turn-off and 
the base of the red giant branch. It has also been suggested that an increase in the surface 
magnetic activity due to evolution can account for the lack of detections in this region, since 
activity depresses the oscillation amplitudes (iGillilandl Il985l ). The concentration of stars 
above the gap (around T^s ~5200 K and i^max ~400 /xHz) further supports the latter idea. 
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Fig. 2. — Position in the log g-Tes plane of the targets, where log g was obtained from 
the scaling relations. Stars with masses determined to be 1-Mq ± 15% are plotted in red 
circles, and stars with oscillations detected that do not match the criterion are plotted as 
gray diamonds (without error bars to reduce clutter). Blue squares are stars with available 
spectroscopic results. Stellar tracks and position of the Sun are plotted as in Fig. [H 
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Fig. 3. — The first three panels show the T^s difference (IRFM- model) for the four targets 
marked in Fig. |2] and: a) the 1-Mq track at solar metallicity, b) a I-Mq track at the 
spectroscopic metallicity for that star, c) a track of the mass determined by the direct 
method and the spectroscopic metallicity for that star. Panel d) shows the log g values 
obtained from spectroscopy compared to those determined with the direct method , with 
the dashed line indicating the one-to-one relation. 
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Fig. 4. — Metallicity distribution in the Kepler field as simulated by the SMAUG code, 
showing stars with a probabihty higher than 90% to have oscillations detected after one 
month of observations. 
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Star 


Mass (Mq)^ 


Teff IRFM (K) 


Teff^ (K) 


log 5^ 


log 


[Fe/H]'^ 


Source 


1 


0.96 ±0.051 


5671 ± 94 


5715 ±82 


4.32 ± 0.007 


4.31 ±0.14 


0.30 ±0.06 


ARES 


2 


1.04 ± 0.043 


6073 ± 100 


6073 ± 78 


4.28 ± 0.006 


4.38 ±0.12 


-0.10 ±0.06 


ARES 


3 


1.07 ±0.051 


6105 ± 101 


5940 ± 70 


4.18 ±0.006 


4.21 ±0.08 


-0.10 ±0.07 


VWA 


4 


1.07 ±0.062 


5836 ± 124 


5870 ± 70 


4.09 ± 0.008 


4.07 ± 0.08 


-0.01 ±0.07 


VWA 



Table 1: Stellar parameters of the four targets with available spectroscopic observations. 



"From the direct method 
''From spectra 



4. Results 



We have applied the direct method to the sample in Fig. [T] and determined log g values 
using the asteroseismic results extracted by the QML pipeline, which has been used a s 
reference for comparison with other pipelines ( IVerner &: Roxburghll2010l . lVerner et al.ll201ll ). 
To identify stars with masses similar to the Sun for constructing the evolutionary sequence, 
we selected targets that according to their determined masses and uncertainty, propagated 
from Eq. [H had a probability larger than 68% to be in the 0.85-1.15 Mq range. 

Figure [2] shows these stars in the log g-T^s plane, where 72 targets are found to match 
the selection criterion. The distribution of matches in the observational plane is neatly 
contained by the stellar tracks, pointing towards an agreement between the asteroseismic 
mass values and evolutionary predictions. As in Fig. [H a gap can be seen around Teg ~5300 
K and log(g)~3.7, again identified with the separation between the end of the main-sequence 
and the red giant phases. 

In Fig. [2] it is also clear that most of the matches fall on the hotter side of the solar- 
metallicity 1 track. There are two possible reasons for this. The first is that the 



am plitudes of the osci l 



see 



Samadi et al. 



2011 



2007, 



ations increase wi t h increasing L IM ratio (or some variant of it 



Mosser et al.ll2010l . IChaplin et al.ll2011bl . iKieldsen fc Beddindll995 



Huber et al.ll201ll ). resulting in a detection bias in our ensemble that favors hotter 
(higher-than-solar-mass) main-sequence and subgiant stars. However, for our considered 
threshold, only 53 % of the matches have a determined mass larger than l-M©, suggesting a 
mass distribution barely skewed to higher values than our target mass. The second reason 
is a spread due to different chemical compositions of the stars because, for a given mass, a 
lower metallicity will imply a higher luminosity (and vice versa). 

Keeping this in mind, we can infer from the position of the stellar tracks in Fig. [2] 
that the spread is most likely due to differences in the chemical compositions of the stars. 
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The effects of different masses and other physical processes such as microscopic diffusion 
(estimated by the separation between the I-Mq solar metallicity track and the position of 
the Sun) are not large enough to account for this spread. Unfortunately, we do not yet 
have spectroscopic determinations of metallicity for most of the stars we are studying, so 
the question arises: can we reliably identify stars with similar masses, or does the direct 
method provide incorrect mass determinations that are camouflaged by the unknown stellar 
parameters? 

We can have a first glimpse at the answer by using spectroscopic results which are 
available for four of our target s (blue squares in Fig. [2]). The stellar pararneters have been 
determined using the ARES JSousa et al.l boosh and VWA JSruntt et al.l boioh methods, 



where we selected the spectroscopic set of parameters with log (7 value closer to the direct 
method determination. These are shown in Table [H It can be seen that our Tgg input values 
(IRFM) for these four stars are compatible with the spectroscopic results. 

In Fig. [3]we compare the results of the direct method for the four stars with evolutionary 
predictions by showing deviations in T^s between our IRFM input parameter and stellar 
tracks at the log g value determined from the direct method. Figure [3^ shows the T^s 
difference between our targets and a I-Mq evolutionary track at solar metallicity, with large 
deviations arising from a combination of different metallicities and masses. The following 
panels show the results when we take these differences into account step by step. Firstly, 
Fig. [3)3 depicts the difference in Tg^ between our targets and 1 Mq tracks at the metallicity 
obtained from spectroscopy. Stars 1 and 2 are already compatible with the stellar tracks, 
which is encouraging as their masses determined with the direct method are consistent with 
I-M0 within their 1-a uncertainties. The other two stars have slightly higher-than-solar mass 
determinations. In Fig. ^ the difference is calculated with respect to evolutionary tracks 
of the metallicity measured by spectroscopy and mass determined from the direct method. 
It is clear that the agreement improves as we include the effects of mass and chemical 
composition, reinforcing the conclusion that the direct method predicts masses for field stars 
compatible with evolutionary tracks of the metallicities obtained from spectroscopy, in a 
model-independent way. Moreover, the log g values determined with the direct method are 
in agreement with the spectroscopic results, as shown in Fig. [311. 

To obtain some information on the expected metallicities for the rest of our targets, 
we simulated the stellar population in the K epler field using the Simple Model for Analytic 
Understanding of our Galaxy code rSMAUG. ISchonrich fc Binnevll2009al lbll The outcome of 



the simulation has been filtered to retain only synthetic stars that, according to their radius, 
effective temperature and visual magnitude, have a probability gre ater than 90% for th e 



Kepler satellite to detect oscillations in one-month-long observations ( IChaplin et al.ll2011bl ). 
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The resulting subset has been further restricted to stars with Tcs < 7000 K, log g > 3.4 and 
rriv < 12.5. The metallicity distribution of that population is shown in Fig. |4] and peaks 
at a sub-solar value, with the bulk of stars contained in a range that makes them clearly 
compatible with the position of our matches and stellar tracks in Fig. |2J 



Discussion 



The sequence of field stars identified in the previous section shows very good agreement 
with the expected position of 1-Mq stars in the log g-T^s plane. For those targets with 
metallicity measurements available, the masses and radii determined with the direct method 
are fully compatible with evolutionary tracks and spectroscopic log g determinations. Never- 
theless, we must keep in mind that uncertainties in the observables and the scaling relations 
can affect our derived quantities. 

From Eqs. [T] and [2] it is clear that the mass and radius determinations are mostly 
sensitive to the uncertainties arising from the asteroseismic parameters. Although the errors 
in effective temperature determinations cannot be neglected, the fractional uncertainties 
in i/max are currently up to four times larger than those in Tes and Au. Therefore, the 
uncertainties in the frequency of maximum oscillation power currently account for most of 
the error budget. 

Masses and radii determined from the direct method have an intrinsic uncertainty asso- 
ciated with two effects: deviations from the scaling relations and unknown stellar properties 
that can change the pulsation characteristics of the star. The scaling relation between the 
large frequency sepa ration and the mean stellar density has been shown by theoretical models 
to be quite robust (IStello et al.ll2009al IWhite et al.ll201l[ ). It is more complicated to esti- 
mate the uncertainties arising from the z/max scaling, since it relies on the relation between 
z/jnax and the acoustic c ut-off frequency under the assumption of an isothermal atmosphere 
(e.g. iBrown et al.lll991[). Although a partia l theoretical basis for this relation has only re- 



cently been developed ( iBelkacem et al. 



201ll ). comparisons under this assumption with stellar 



models predict tha t it holds within a few percent for cool models (close to solar-mass, see 
Stello et al.ll2009al ). Observational results for stars whose parameters are accurately known 
confirm this result ( iBedding fc KjeldsenI 120031 ). Realistic model atmospheres could help to 
assess its deviations, as could comparisons with direct observational methods to measure 
stellar properties, such as long-baseline interferometry. 



The other main source of uncertainties is related to unknown stellar parameters, such 
as chemical composition. It is possible to complement the direct method with our knowledge 
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of stellar evolution and match the seismic observables to their expected values from previ- 
ously calculated evolutionary models. This so-called grid-based method has be en shown by 
simula t ions to perform as well as the direct me thod for constraining the radius ( IStello et al. 



2009bl . iBasu et al.l l2010l . iKallinger et al.ll2010[ ). When accurate metallicity measurements 



are available, th e grid-based method can reduce the uncertainty in the mass determination 
( iGai et al.ll201ll ). A combination of both methods will allow us to increase the precision and 
make the results even more robust. Independent t ests are currently bein g developed to assess 
their uncertainties (see, for example, Sect. 3.3 in IChaplin et al.ll2011al ). Another validation 
of our results can be made by comparing the log g values determined with the direct method 
to the parallaxes that will be obtained from the Kepler and Gaia missions in the future, 
give n that a variation in mas s of 10 % translate into a change in logg of approximately 0.04 
dex ( ICasagrande et al.l 1201 if ). 



In this work, we constructed the evolutionary sequence by applying the selection cri- 
terion to the parameters extracted by one particular pipeline. Different pipelines produced 
similar results in terms of the distribution of the matches in the HR diagram and the spread 
around the evolutionary tracks (see Fig. [2]). A natural step forward would be to combine the 
results of several pipelines to restrict further the matches and make the results even more 
robust. However, there is still some variation between the extraction methods in the deter- 
mination of z/jnax- The reason probably lies in the low signal-to-noise ratio o f some of the 
obser vations and the relatively short time-span of the data currently available (jVerner et al. 
201lh . 



Finally, our selection criterion can be arbitrarily modified to accept, for instance, stars 
with masses within 10% of 1.0 M0. In such case only 20 targets satisfy the criterion, while 
increasing the threshold to 20 % results in 120 matching stars. We have kept the 15 % 
threshold as a good compromise between a required high precision and the uncertainties 
described in the previous paragraphs. 



6. Conclusions and further development 

Constructing empirical evolutionary sequences from asteroseismic data offers the excit- 
ing possibility of performing differential analysis on field stars of similar masses, allowing us 
to test stellar properties very precisely through different evolutionary phases. Scaling rela- 
tions using global asteroseismic parameters and T^s measurements allowed us to determine 
masses of stars in a model-independent way. The position of our matches in Fig. [2] suggests 
that we have successfully identified an evolutionary sequence of field stars with masses very 
close to 1 Mq. The results are encouraging and for the first time we can construct such a 
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sequence without parallax information, spectroscopic log g, or masses estimated from evo- 
lutionary tracks. Longer time series and ground-based follow-up spectroscopy will further 
enhance the capabilities of this technique, allowing us to reach a higher level of precision. 

The results from Sect. H] indicate that our mass and radius determinations are very ro- 
bust. Using these results and the input effective temperature, together with a grid of stellar 
models, it should be possible to derive the expected metallicity distribution of the targets. 
This distribution will depend on the input parameters used to calculate the grid, most no- 
tably the selected solar abundances, AY/AZ value, aMLT parameter, and the inclusion (or 
not) of mixing processes such as diffusion. Future comparisons of this metallicity distribu- 
tion with accurate chemical composition determinations could help to constrain the input 
parameters of the grid models that produce the largest variation in the resulting metallicities. 

When suitable analogues are found in different evolutionary stages and their individual 
frequencies determined, a differential analysis can be performed to better constrain the inte- 
rior physics of stars throughout their evolution. If targets with masses higher than solar are 
chosen instead, we could differentially study the effects of mixing processes in their interiors 
and constrain processes that still rely mostly on empirical calibrations, such as convectiv e 



core evolution and overshooting (e.g. iDe Meulenaer et al.l l2010l ISilva Aguirre et al.ll201ll ). 



For example, by selecting a sequence of stars of very similar mass it could be possible to 
produce a precise relative age calibration. 

Asteroseismology has provided us with a method to disentangle stars still ascendin g 



the red giant branch from those already burning helium in their cores ( iBedding et al.ll201ll ). 
Our work can be extended to identify evolutionary sequences from the main-sequence to 
the red clump. Differential analysis of these targets can provide further constraints on the 
main-sequence physical processes and shed light on the progenitors of horizontal branch 
stars. 

In the near future, it will be interesting to construct evolutionary sequences and select 
twins of very well-studied main-sequence stars, both seismically and in terms of their stellar 
parameters. From the ensemble it is already possible to identify sequences for a Cen A and B, 
although for the latter case the number of analogues is much smaller due to our bias towards 
higher amplitudes than the Sun. Some twins of a Cen A are already present in the ensemble, 
and longer time-series will allow us to expand the search to lower temperatures and find twins 
of our Sun and a Cen B. The potential of combining asteroseismic dat a with spectroscop ic 



determinations for solar twins has only recently begun to be exploited (IBazot et al.ll201ll ). 
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